Magnetic moment of the A(1232)-resonance in chiral effective field theory 
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We perform a relativistic chiral effective field theory calculation of the radiative pion photopro- 
duction (jp — + 7r''p7') in the A-resonance region, to next-to-leading order in the "delta-expansion". 
This work is aimed at a model-independent extraction of the magnetic moment from new precise 
measurements of this reaction. It also predicts the chiral behavior of A's magnetic moment, which 
can be used to extrapolate the recent lattice QCD results to the physical point. 

PACS numbers: 12.39.Fe, 13.40.Em, 25.20.Dc 



The A(1232)-isobar is the niost distinguished and weU- 
studied nucleon resonance. However, such a fundamen- 
tal property as its magnetic dipole moment (MDM) has 
thusfar escaped a precise determination. The problem 
is generic to any unstable particle whose lifetime is too 
short for its MDM to be measurable in the usual way 
through spin precession experiments. A measurement of 
the MDM of such an unstable particle can apparently be 
done only indirectly, in a three-step process, where the 
particle is first produced, then emits a low-energy photon 
which plays the role of an external magnetic field, and fi- 
nally decays. In this way the MDM of A++ is accessed 
in the reaction n'^p — > ir'^pj 0,0 while the MDM of A+ 
can be determined using the radiative pion photoproduc- 
tion (7p TT^jry') 

A first experiment devoted to the MDM of A+ was 
completed in 2002 The value extracted in this ex- 
periment, /LtA+ = 2.7^i['3(stat.) ± 1.5(syst.) ± 3(theor.) 
[nuclear magnetons], is based on theoretical input from 
the phenomenological model H, 0] of the jp — > n^pj' 
reaction. To improve upon the precision of this mea- 
surement, a dedicated series of experiments has recently 
been carried out by the Crystal Ball Collaboration at 
MAMI 0. These experiments achieve about two or- 
ders of magnitude better statistics than the pioneering 
experiment The aim of the present work is to comple- 
ment these high-precision measurements with an accu- 
rate and model-independent analysis of the jp ir^p'y' 
reaction, within the framework of chiral effective field 
theory (xEFT). 

The xEFT of the strong interaction is indispensable, 
at least at present, in relating the low-energy observ- 
ables (e.g., hadron masses, magnetic moments, scatter- 
ing lengths) to ab initio QCD calculations on the lattice. 
On the other hand, xEFT can and should be used in 
extracting various hadronic properties from the experi- 
ment. In this Letter we will show how xEFT fulfills both 
of these roles in a gratifying fashion. The one-loop cal- 
culation we present here is sufficient to both complete 



the next-to-leading order calculation of the 7p — > ir'^pj' 
reaction, in the A-resonance region, and perform a chiral 
extrapolation of lattice QCD results for A's MDM [1,0]. 

Our starting effective Lagrangian is that of the chi- 
ral perturbation theory (^PT) with pion and nucleon 
fields 01 . The A then is included explicitly in the so- 
called (5-expansion scheme 0|. We organize the La- 
grangian £^*\ such that superscript i stands for the 
power of electromagnetic coupling e plus the number of 
derivatives of pion and photon fields. Writing here only 
the terms involving the spin-3/2 isospin-3/2 field of 
the A-isobar gives: ^ 
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where AI and Ma are, respectively, the nucleon and 
A-isobar masses, N and tt°- (a = 1,2,3) stand for the 
nucleon and pion fields, Z?^ is the covariant derivative 
ensuring the electromagnetic gauge-invariance, F'^'^ and 
F'^'^ are the electromagnetic field strength and its dual, 
Ta are the isospin 1/2 to 3/2 transition matrices, cou- 
plingconstants U = 92.4 MeV, gM = 2.94, gs = -0.96, 
see [nl [13 for further details. The MDM is defined 



^ Here we introduce totally antisymmetric products of 7-matrices: 
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here in units of [e/2M/s\- We omit the higher electro- 
magnetic moments, because they do not contribute at 
the orders that we consider. 

Note that £a2, contains the free Lagrangian, which is 
formulated in |l3| such that the number of spin degrees 
of freedom of the relativistic spin-3/2 field is constrained 
to the physical number: 2s + 1 = 4. The to A tran- 
sition couplings in Eq. are consistent with these con- 
straints [ij, |l3 • The 7A A coupling is more subtle since 
in this case constraints do not hold for sufhciently strong 
electromagnetic fields, see, e.g., We do not deal 

with this problem here, thus assuming the electromag- 
netic field to be weak, compared to the A mass scale. 

We now briefly describe the power counting in the 
i5-expansion scheme. The excitation energy of the A- 
resonance, i.e., A = A/a ~ Mm — 293 MeV is treated 
as a light scale, so A <c A, where A ~ 1 GeV stands 
for the heavy scales of the theory. At the same time, 
A is counted differently from the other light scale of the 
theory - the pion mass, m^r. Namely, A/ A counts as 
one power of the small parameter 5 while m-^j h. counts 
as two powers of 8. Each graph can then be character- 
ized by an overall i5-counting index n, which simply tells 
us that the graph is of size 5". Because the theory has 
two distinct light scales (mjr and A) the i5-counting in- 
dex depends on whether the characteristic momentum 
p is in the low-energy region (p m^) or in the res- 
onance region (p ~ A). In the low-energy region the 
index of a graph with L loops, TV^ pion propagators, 
'Njsi nucleon propagators, A^a A-isobar propagators, and 
Vi vertices of dimension z is n = ^n^^T — N/^, where 
n^PT = J2t iL~ Nn ~ 2N^ is the index in xPT 

with no A's [ig. 

In the resonance region, one needs to distinguish the 
one-A-reducible (OAR) graphs, because they contain A 
propagators which go as l/(p— A) and hence such graphs 
are large and all need to be included. Their resummation 
amounts to dressing the A propagators so that they be- 
have as l/{p— A — S). The self-energy S begins at order 

and thus, for p ^ A, the dressed OAR propagator goes 
as If the number of such propagators in a graph 

is Nqar, the power-counting index of this graph in the 
resonance region is given by n = n^PT — A^A — ^Nqar- 

Consider now the amplitude for the jp — *■ ir^p'y' re- 
action. The optimal sensitivity to the MDM term is 
achieved when the incident photon energy is in the vicin- 
ity of A, while the outgoing photon energy is of order of 
m,r- In this case the 7p — > n^pj' amplitude to next-to- 
leading order (NLO) in the 5-expansion is given by the 
diagrams ^a), (b), and (c), where the shaded blobs, in 
addition to vertices from Eq. contain the one- loop 
corrections shown in Fig. ^d), (e), (f). To present our 
calculation of the loops we introduce the mass ratios: 
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FIG. 1: Diagrams for the 7p iv^p-f' reaction at NLO in the 
(5-expansion, considered in this work. Double lines represent 
the A propagators. 

/Ma, r = M/Ma. The self-energy, Fig.Utd), 
has the following form 

S^^'i^) = A{p^)j^'"^p^ + B(p2)7''^ (2a) 

where the scalar functions, after dimensional regulariza- 
tion, take the form 

A(p^) = -^C^ f dxxM"^ {L-l + lnM"^), (2b) 

S(p2) -\C^r ( dxM'^ {L-l + \nM^), (2c) 
"'0 

with C = hAMA/iSnU), L = -2/(4 - rf) + 7s + 
ln(47rAfA/A), c? ^ 4 the number of dimensions, je = 
— F'(l) ~ 0.5772, A the renormalization scale and 

M'^ix) = xy?' + xy - x{l - x){p^/Ml) - ie. (3) 

After the on-mass-shell wave-function and mass renor- 
malizations the NLO A-propagator is given by 

'""^^^ (^- Ma)[1 - ilmS'iMA)] - ilmS{MA) ' 

(4) 

where we have introduced the spin-3/2 projec- 
tion operator, 7-'(3/2)^ ^^^^^ ^^xe following definitions: 
S{Ma) = MaA{MI) + B{Ml), S'{Ma) = A{Ml) + 
2MA{d I dp'^)[MAA{p^) + B{p^)]p2^Ml- These functions 
are complex when the A-isobar is heavier than the pion 
production threshold. A/a > M + m^r. In this case 
A^^ < and hence the logarithm in Eq. jSJ gives rise 
to an imaginary part: 

Imi:(MA) = -(27r/3)MAC2(cH-r) A^ (5a) 
Imi;'(MA) = -2ttC'^ \[a{l - a){a + r) 

- \\\r + r^-^l'')\, (5b) 

with a = i(l + r^ — /i^), A = — r^. We note that the 
width of the resonance is given by Fa = — 2Imi7(AfA), 
and find that the experimental value Fa — 115 MeV 
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translates into Ha — 2.85, the value which we shall use 
in the numerical calculations. 

The 7AA vertex, omitting the electric quadrupole and 
magnetic octupole terms, is written in the form: 



= eua{p' 



where e'' is the photon polarization vector, u" is the 
vector-spinor of the A. In this notation, the MDM is 
given by /iA = -P(O). The Ward-Takahashi identity, 

q^r^''^'{p',p) - e [{S-'r^'{p') - {S-'r^{p)] , (7) 

demands that F{0) + G{0) = 1-I:'{Ma). This condition 
is verified explicitly in our calculation. Since we have al- 
ready given the expression for the self-energy, we present 
here only the expressions for G. The contributions of 
diagrams 1(e) and at = 0, are: 

G('='(0) = --C^ [ dxx{l - 2x){x - r) 
Jo 

x{L + ln[xfi^ + (1 - x)r^ - x{l - x) - ie]} , (8) 

G^^) (0) = -2^2 /" dxx^{l-x- r) 
Jo 

X {L + ln[xr^ + (1 - x)fi'^ - x{l ~ x) - ie]} . (9) 

Their contribution to /iA for different isospin states is, 
e.g., fi^llf = F(-'^+F(f\ = il/3)F(') + i2/3)F^f\ 
As all lattice data for fXA at present and in the foresee- 
able future are for larger than the physical values of m,r , 
their comparison with experiment requires the knowledge 
of the m,r-dependence for this quantity. In contrast to 
the heavy-baryon result our xEFT calculation is 

manifestly relativistic, and, as argued earlier for the nu- 
cleon Til, should be better suited for such extrapolations 
of lattice data. Figure[21shows the pion mass dependence 
of real and imaginary parts of the A+ and A++ MDMs, 
according to our one-loop calculation. Each of the two 
solid curves has a free parameter, the counterterm //a 
from C\ , adjusted to agree with the lattice data at larger 
values of m^. As can be seen from Fig. 12 the A MDM 
develops an imaginary part when < A = Ma — M, 
whereas the real part has a pronounced cusp at m^r = A. 
For /iA+ our curve is in disagreement with the trend 
of the recent lattice data, which possibly is due to the 
"quenching" in the lattice calculations. The dotted line 
in Fig. 121 shows the result |0| for the magnetic moment 
for the proton. One sees that /xa+ and fXp, while having 
very distinct behavior for small 711,^, are approximately 
equal for larger values of TOtt. 
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FIG. 2: Pion mass dependence of the real (solid curves) and 
imaginary (dashed curves) parts of A^""" and A""" MDMs [in 
nuclear magnetons]. Dotted curve is the result for the pro- 
ton magnetic moment from Ref. [T^ . The experimental data 
point for A"'"*' is from PDG analysis [l9l| . Lattice data are 
from H for A++ and from 91 for A+. 



We next discuss our results for the jp — > ir^pr)' ob- 
servables. The NLO calculation of this process in the 
(5-expansion corresponds with the diagrams of Fig. ^ 
As outlined above, this calculation completely fixes the 
imaginary part of the 7AA vertex. It leaves /iA as only 
free parameter, which enters as a low energy constant in 
L'^'^\ Thus the real part of //a+ is to be extracted from 
the 7p — *■ n'^p'y' observables, some of which are shown in 
Fig. 121 for an incoming photon energy = 400 MeV 
as function of the emitted photon energy E'^^-"^- . In the 
soft-photon limit {E'^'^-"^- 0), the 7p ir^p^' reaction 
is completely determined from the bremsstrahlung pro- 
cess from the initial and final protons. The deviations of 
the 7p ~> tt'^pY observables, away from the soft-photon 
limit, will then allow to study the sensitivity to /iA+- It 
is therefore very useful to introduce the ratio 0: 



R = 



1 



E' 



da 



''dEL 



(10) 



where du/dE'^ is the 7p n'^pj' cross section integrated 
over the pion and photon angles, and a^^ is the angular in- 
tegrated cross section for the jp — *■ n'^p process weighted 
with the bremsstrahlung factor, as detailed in This 
ratio R has the property that in the soft-photon limit, the 
low energy theorem predicts that i? — > 1. From Fig. |3 
one then sees that the EFT calculation obeys this theo- 
rem. This is a consequence of gauge-invariance which is 
maintained exactly throughout our calculation, also away 
from the soft-photon limit. 

The EFT result for R shows clear deviations from unity 
at higher outgoing photon energies, in good agreement 



4 



E.!"'' = 400 MeV 




0.8 - 
0.6 i 



e'"'' = 400 MeV 

e/'" =90° 



0.4 - 

: Re\l^ = 5 

0-2 : Rep.l = l 

Re = 3, Im (i^ = 

\ 

0.2 - o°<e;"' <90° 
: 0° < <5^"" < iso° 

.y 0.1 - 




20 40 60 80 100 120 140 160 



E„ (MeV) 



FIG. 3: The outgoing photon energy dependence of the 
7P Tv'''p'y' observables for different values of (in units 
e/2AfA). Top panel: the ratio of 7p 7v"pj' to jp — > yr^p 
cross-sections Eq. Hl()|l . Data points are from 4,]. Mid- 
dle panel: the linear-polarization photon asymmetry of the 
7P Tv'^p-y' cross-sections differential w.r.t. the outgoing pho- 
ton energy and pion cm. angle. The data point ai E!y — 
corresponds with the jp yr^p photon asymmetry from f20| . 
Lower panel: the circular-polarization photon asymmetry (as 
defined in 0), where the outgoing photon angles have been 
integrated over the indicated range. 



with the first data for this process The sensitivity 
of the EFT calculation to the /xa is a very promising 
setting for the dedicated second-generation experiment 
which has recently been completed by the Crystal Ball 
Coll. at MAMI [t'I. It improves upon the statistics of the 
first experiment (Fig. ^ by at least two orders of magni- 
tude and will allow for a reliable extraction of fi^+ using 
the EFT calculation presented here. 

Besides the cross section, the — s- ir'^pj' asymme- 
tries for linearly and circularly polarized incident pho- 
tons have also been measured in the recent dedicated 
experiment 0|- They are also shown in Fig. The 
photon asymmetry for linearly polarized photons, S, at 
E'^ = exactly reduces to the 7p ir^p asymmetry. It 



is seen from Fig. |31 that our calculation is in good agree- 
ment with the experimental value. At higher outgoing 
photon energies, the photon asymmetry as predicted by 
the NLO EFT calculation remains nearly constant and is 
very weakly dependent on /xa- It is an ideal observable 
for a consistency check of the EFT calculation and to 
test that the A diagrams of Fig. ^ indeed dominate the 
process. Mechanisms involving 7r-photoproduction Born 
terms followed by TriV rescattering have been considered 
in model calculations 0, @ . In the 5-counting they start 
contributing at next-next-to-leading order and therefore 
will provide the main source of corrections to the present 
NLO results. 

The asymmetry for circularly polarized photons, Scire > 
(which is exactly zero for a two body process due to 
refiection symmetry w.r.t. the reaction plane) has been 
proposed [6j as a unique observable to enhance the sen- 
sitivity to /ZA- Indeed, in the soft-photon limit , where 
the 7p TT^pj' process reduces to a two-body process, 
Ticirc is exactly zero. Therefore, its value at higher outgo- 
ing photon energies is directly proportional to /iA- One 
sees from Fig. |31 (lower panel) that our EFT calculation 
supports this observation, and shows sizeably different 
asymmetries for different values of /ua • A combined fit of 
all three observables shown in Fig. O will therefore allow 
for a very stringent test of the EFT calculation, which 
can then be used to extract the A+ MDM. 



In conclusion, we have performed a manifestly gauge- 
and Lorentz-invariant chiral EFT calculation for the 
7P —> tt'^p^' reaction in the A(1232) resonance region. 
To next-to-leading order in the (5-expansion, the only 
free parameter entering the calculation is the A+ mag- 
netic dipole moment /xa+- Due to the unstable nature 
of the A-isobar its magnetic moment acquires an imagi- 
nary part, an effect which is computed in this work and 
will be taken into account in the extraction of /iA from 
experiment. Our present calculation is found to be in 
good agreement with first experimental results for the 
7P TT^p'y' cross sections, and will allow for a model- 
independent extraction of /LtA+ from a combined fit to 
cross sections and photon asymmetries measured in new 
dedicated experiments. At the same time, this chiral 
EFT calculation provides a crucial connection of present 
lattice QCD results for /ia at values of m^r > Ma — M 
to the physical pion mass. 
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